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Using gauge/gravity duality, we study emission and absorption rates of scalar quanta from far- 
from-equilibrium N = 4 supersymmetric Yang-Mills plasma. By comparing the emission and absorp- 
tion rates to expectations from the Fluctuation-Dissipation Theorem, we study how the spectrum 
thermalizes and how the thermalization time depends on the four momentum of the emitted quanta. 



Introduction-Gauge/ gravity duality [1], or hologra- 
phy, is a powerful tool for studying real-time dynamics in 
strongly coupled quantum field theories. Through holog- 
raphy the creation and thermalization of non-abelian 
plasma maps onto the process of gravitational collapse 
and black hole thermalization, which can be studied nu- 
merically. The equilibration of holographic plasma can 
yield valuable insight into the thermalization of Quark- 
Gluon Plasma (QGP) produced in heavy ion collisions. 
In addition, precise questions about thermalization in the 
quantum field theory can clarify the statistical properties 
of non-equilibrium black holes. The simplest theory to 
study with a gravitational dual is AT = 4 supersymmet- 
ric Yang-Mills theory (SYM) at large N c and 't Hooft 
coupling A. 

As analog problem of photon production in the QGP, 
we study the emission of weakly coupled (Ad) dilatons 
from strongly coupled non-equilibrium SYM plasma. 
Dilaton emission and absorption rates are proportional 
to the Fourier transform of the Wightman functions, 
9-{x\x') = (O(x')d(x)) and g+(x\x') = (6(x)d(x')) re- 
spectively, where 0(x) is the SYM Lagrange density. In 
thermal equilibrium at temperature T these rates are re- 
lated by the Fluctuation Dissipation Theorem (FDT) 

g-(q) = e-»/ T g + (q), (1) 

where q = (u>, q) is the four-momentum conjugate to 
x — x'. The FDT is a direct consequence of the fact that 
the density matrix of the system takes its equilibrium 
form e~ H l T . Therefore, the degree to which the FDT is 
satisfied in non-equilibrium plasmas provides an unequiv- 
ocal measure of the equilibration a particular Fourier 
mode q. 

Through holography, the emission and absorption of 
weakly coupled quanta from non-equilbrium SYM plasma 
is dual to the emission and absorption of Hawking ra- 
diation from a non-equilibrium 5d black brane. More 
precisely, the Ad SYM Wightman functions g± are en- 
coded in 5d Wightman functions G± computed in the 
black brane geometry. The Wightman functions satisfy 
classical equations of motion which can be solved numeri- 
cally in non-equilibrium geometries after the appropriate 
quantum mechanical initial conditions are specified [2]. 
Only the absorption rate of a classical Ad dilaton field, 



g + — g_, is determined by the classical response of the 
black brane [3, 4]. 

For definiteness, we will study the thermalization of 
homogeneous but anisotropic SYM plasma. To create 
this non-equilibrium state, we consider the response of an 
initially equilibrium SYM plasma at temperature 7] to a 
time-dependent gravitational field [5] . The homogeneous 
Ad metric of the external field is 

ds 2 = -dt 2 + e b ^ dx\ + er 2b ^ dx 2 , (2) 

where x± = {0:1,2:2}, and the time dependence of the 
field is Gaussian b(t) = c/(V27rT 2 )e-* 2 / 2r2 . The time- 
dependent geometry does work on the SYM plasma and 
creates a non-equilibrium state with an anisotropic stress 
(rpn"} — diag[£, V±, V±, V\\], which subsequently equli- 
brates at late times. The amplitude and width of b(t), 
c = 2.11, T = 1/irTf, are adjusted using the simulations 
described below so that the energy density changes by a 
factor of fifty, and the final temperature is Tf = 50 1//4 Ti. 

We will quantify and analyze the equilibration of sys- 
tem by monitoring the emission and absorption rates in 
a given frequency band as a function of time. Clearly, 
the frequency and time resolutions of this analysis are 
limited by the uncertainty principle, AuAt > ™. To 
achieve the best possible resolution in both frequency and 
time we will compute Wigner transforms of g±(x, x') av- 
eraged over phase-space with a minimal uncertainty wave 
packet. The Wigner transforms are 

g±(t,e>,q) = J d(Ax)g ± (x\x')e~^ Ax , (3) 

with Ax = x — x' and i = |(i + i 7 ), an d the minimally 
averaged transforms, known as Gabor transforms, are 

«7 ± (*>,g) = / ^ e -^ e -^ fc(? y )?)i 

(4) 

where Auj = uj' — u, At = v — t, and a is a temporal res- 
olution parameter which we choose to be a — l/nTf. For 
the homogeneous system under study, the energy den- 
sity is constant once the external metric is static, and 
an effective temperature can be defined via the energy 
density 

-1- = T cS (t) = {8\£(t)\/(3n 2 N 2 )) 1/4 . (5) 
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If 9± (<?) satisfies the FDT at late times, then the Gabor 
transforms satisfy 

(6) 

The degree to which (6) is satisfied provides an excellent 
definition of equilibration. 

Gravitational description. — The 5d gravitational dual 
to an anisotropic SYM plasma created by a time depen- 
dent Ad geometry was studied in [5] . We follow the anal- 
ysis of [5] and outline some of the salient features. 

Following [5], using translation invariance we may 
write the 5c? metric as 

ds 2 = -Adt 2 + Y? [e B dx\ + e- 2B dxf { ] + 2dr dt , (7) 

where A, B, and £ are functions of r and t, only. The co- 
ordinates t and r are generalized Eddington-Finkelstcin 
coordinates. Infalling radial null geodesies have constant 
values of t (as well as x± and xm). The Ad boundary of 
the 5c? geometry is located at r = oo, and is proportional 
to the geometry of the SYM plasma in eq. (2) . 

The evolution of the 5c? metric (7) is determined by 
solving numerically Einstein's equations with negative 
cosmological constant [6] . Demanding that the boundary 
metric is that of (2), equates to imposing the boundary 
conditions, lim,--^ B(t, r) = b(t), lim r _ s . 00 £(f, r)/r = 1 
[5]. The near-boundary behavior of the metric deter- 
mines the expectation value of the SYM stress (T^ v ) [7]. 

In the infinite past, when b(t) = 0, we take the initial 
geometry to be an equilibrium black brane of tempera- 
ture T;. When the boundary source turns on, b'(t) ^ 0, 
the changing boundary geometry creates gravitational 
radiation that propagates into the bulk creating a non- 
equilibrium black brane. When the boundary geometry 
is static, b'(t) = 0, the 5d geometry relaxes to an equilib- 
rium black brane with temperature T{. After the solving 
the Einstein equations, the event horizon r^it) can be 
found by tracing null radial geodesies backwards in time 
from the far future. 

As the black brane evolves and equilibrates, it contin- 
uously emits and absorbs Hawking radiation. The radi- 
ation propagates up from the horizon to the boundary, 
where it is reflected back towards the horizon and is ab- 
sorbed. The process continues until the 5c? dilaton modes 
are in equilibrium with the black brane and are thermally 
occupied with distribution function n(u) = l/(e"/ Tt — I). 
This competitive evolution between emission and absorp- 
tion is the gravitational dual of thermalization in SYM 
plasmas. 

The emission of 5d dilatons is encoded in the 5c? Wight- 
man function G-{x\, r\\x 2 , r 2 ) = (<p(x 2 ,r 2 )ip(xi,r 1 )) 
where tp is the dilaton field. Similarly, the absorp- 
tion of 5c? dilatons is encoded in G+(xi,rx\x 2 ,r 2 ) = 
(ip(xi, ri)ip(x2, fi))- Both G± satisfy the classical equa- 



tions of motion 

- Df 1) G ± (x 1 ,r 1 \x 2 ,r 2 ) = -D 2 2 - ) G±(xx,r 1 \x2,r 2 ) = 0, 

(8) 

where -Df n ) is the covariant wave operator under 
the background geometry with respect to coordinates 

Previously, we (together with S. Caron-Huot) formu- 
lated a tractable numerical recipe for solving eqs. (8) [2]. 
Our prescription is based on the following logic. First, 
we note that G± satisfies the homogeneous equations of 
motion (8) and is not a Green function. Thus, after spec- 
ifying the initial conditions in infinite past, forward evo- 
lution determines G± in the future. In the equilibrium 
black brane geometry of the distant past, most of the ini- 
tial data specified for G± falls through the black brane 
horizon in a finite time and thus becomes irrelevant for 
future evolution. Examining the geodesies of the equi- 
librium black hole geometry, one concludes that the only 
initial data of relevance to the evolution of G± in the 
future is determined by the coincident point singularities 
of the two point functions arbitrarily close to the event 
horizon, i.e. ultraviolet vacuum fluctuations. These UV 
vacuum fluctuations propagate along horizon skimming 
geodesies and eventually emerge (redshifted) from the 
stretched horizon at = rh(t) + e. The flux of data 
passing through the stretched horizon can be traded for 
effective sources defined on the stretched horizon and the 
solution to (8) can be written in terms of the convolution 
integral 

G±{x 1 ,r 1 \x 2 ,r 2 ) = / {g k (x x , ri\x[, r[)G R (x 2 , r 2 \x 2 , r 2 ) 

x^(z'il4)}, (9) 

where the horizon correlators G± are the effective sources 
and the retarded Green's function Gr satisfies 

n2 ^ / i \ S 4 (x 1 -x 2 )S(r 1 -r 2 ) 

- D, 1) G R {x 1 ,r 1 \x 2 ,r 2 ) = , (10) 

V -9{xi,ri) 

with g the determinant of the metric. 

The horizon correlators can be computed analytically 
[2] . Specifically, the propagation of the UV vacuum fluc- 
tuations from arbitrarily close to the event horizon to the 
stretched horizon can be treated using geometric optics. 
Furthermore, the UV vacuum fluctuations generated by 
coincident point singularities are universal and indepen- 
dent of initial conditions for G± . Following the procedure 
of [2] , in the e — > limit we find 

[g(t,r h )g(t',r h )}^ 4 \ 

-^Mt>sch 2 ^^l, (11) 
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where = \d r A{t,r)\ r=Th , r(i) = f dt'K,(t'), and 

-ftT is the normalization of the dilaton action, S — 
~ T / \/~9(® ( l } ) 2 ■ Once the bulk correlators are deter- 
mined, the corresponding SYM correlators are given by 
[2] 



g±(x\x') 



lim 16rV 4 G ± (x,r|a;',r'). 

r,r 7 — >oo 



(12) 



The numerical calculation of G± and hence g± can be 
greatly simplified by noting that the horizon correlators 
can be written 



Ql{x\x') = K 



K(t)K(t')5 3 {x-x') 

2vdv 



2tt 



n ± (2TT^)e- tv{T(t ^ t '^, (13) 



where n_(x) = l/(exp(x) — 1) and ra + = 1 + ?i_ are 
the Bose occupation and stimulation factors respectively. 
From eqs. (9) and (10) it follows that the bulk correlators 
are 

' Ivdv 



G ± (x,r\x',r') 



K 



2tt 



-n ± {2-Kv)F v (x, r)FZ(x',r'), 

(14) 



where F v (x,r) satisfies 



D 2 F u (x,r) 



n{t)e- tl 



r(t) 



-fl(*,r h )]V* 



<5(r-n)<5 3 (a;-a;'). (15) 



Using translation invariance in the spatial directions, 
we introduce a spatial Fourier transform and solve 
eq. (15) for the mode functions F v (t,q,r). The required 
boundary conditions needed to solve eq. (15) are that 
F v (t, q,r) —> at the boundary and that F v {t, q, r) is reg- 
ular at the horizon. The initial conditions for F v (t,q,r) 
in the infinite past can be computed from an equilibrium 
analysis of eq. (15) for a static black brane at temper- 
ature T;. Finally, we note that using field redefinitions 
it is possible to take the — > limit in our numeri- 
cal analysis. Further numerical details will be presented 
elsewhere. 

Discussion. — Fig. 1(a) shows a plot of the SYM stress 
T^ v = diag(— £ ,V±,V±,V\\) as a function of time. In 
the distant past the stress is static and corresponds to a 
low temperature equilibrium plasma of temperature T;. 
When the boundary geometry (2) starts to change, work 
is done on the system, the energy density generally grows 
and the pressures oscillate. In Fig. 1 the shaded region 
indicates that the energy density is still changing, and 
work is being done on the system by the external bound- 
ary source. After the energy density reaches its final 
value, the transverse and longitudinal pressures V± and 
V\\ equilibrate near time, t ~ 3/irTf. 

Fig. 1 (b) shows the emission and absorption rates from 
non-equilibrium plasmas for timelike and lightlike dila- 
tons with w = 8ttT{. More precisely, the solid lines indi- 
cate Gabor transforms of the emission rate g~- (eqs. (3) 
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FIG. 1: (a) The SYM stress tensor T M " relative to the fi- 
nal energy density £f as a function of time. The shaded 
band indicates when the energy density is changing due to 
the work done by the external gravitational field in the gauge 
theory, (b) A non-equilibrium emission rate, cj-(i,q)/g { l nsl (q) 
(see eq. (4)). The emission rate is exhibited for time-like mo- 
menta, with u) — 8nTf and q = 0, and for light-like momenta, 
with q± — §|| = w/v2 and cj = 87rTf. The dashed lines show 
the FDT-expectation for the emission rate, i.e. the rate de- 
rived using the absorption rate p+ (t, q) and the FDT result 
(eq. (6)). (c) The emission rate for time-like momenta (as in 
(b)) and various values of the frequency, (d) The emission 
rate for light-like momenta (as in (b)) at various values of the 
frequency. 
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and (4)), while the dashed lines indicate an estimate for 
the emission rate, g™ T , which uses the absorption rate 
g + and the coarse-grained FDT (eq. (6)). From Fig. 1(b) 
we see that that the emission rates equilibrate only after 
the stress tensor and the absorption rates have equili- 
brated. 

The frequency dependence of the emission rates is 
studied in Fig. 1(c) and (d). The timelike emission rates 
in Fig. 1 (c) approach a universal curve at large frequency, 
indicating that the emission of timelike dilatons equili- 
brates in a finite time as to — > oo. By contrast, the light- 
like emission rates in Fig. 1(d) take an increasingly long 
time to equilibrate as uj — > oo. 

The qualitative features of the equilibration process 
can be understood by examining the physics of non- 
equilibrium Hawking radiation recorded in eq. (9). The 
Hawking flux from the black hole membrane can be con- 
sidered equilibrated when the horizon sources recorded 
by Q\ satisfy the horizon FDT [2-4]. In general, mem- 
brane equilibrium (as defined by the FDT) is reached 
after the background geometry and boundary stress 
have stopped changing. At this point, an equilibrated 
Hawking flux is transmitted to the boundary through 
the retarded horizon-to-boundary Green functions in 
eq. (9). Consequently, the equilibration of the bound- 
ary two point functions is causally delayed relative to 
the boundary stress by the horizon-to-boundary transit 
time, At ~ 2/irT. This explains the slight time shift 
between Fig. 1(a) and Fig. 1(b). 

Similarly, the different equilibration patterns of time- 
like and lightlike modes can be understood with the dual 
physics of Hawking radiation. The emission from the 
black hole membrane is the same for lightlike and time- 
like modes, since eq. (11) is independent of q. Thus, 
the timelike-spacelike difference in the boundary the- 
ory is due to the horizon-to-boundary Green functions. 
For large frequencies, the salient features of these Green 
functions can be understood with geometric optics. In 
addition, for coordinate time t > 2.57rTf the metric is 
relatively close to its equilibrium form, and equilibrium 
geodesies provide qualitative insight to the geometric op- 
tics of the full geometry. Fig. 2 shows a sample of 5d equi- 
librium geodesies that describe the transport of Hawking 
radiation from the stretched horizon to the boundary for 
a specified \q\/u) in the boundary theory [8, 9]. 

Examining Fig. 2, we see that the Hawking radia- 
tion which first arrives at the boundary is character- 
ized by high frequency, uj — > oo with q fixed. These 
modes are transported along null radial geodesies with 
a minimal transit time of order At ~ 2/wT, and thus 
we expect timelike modes to thermalize first. Similar 
conclusions were also reached in [10]. By contrast, the 
geodesies that determine the Green functions of hard 
lightlike modes have uj ~ \q\ large, and small virtu- 
ally (w 2 — \q\ 2 )/ui 2 ~ \/(uia), where a = 1/ttT { is 
the smearing width in eq. (4). These geodesies (which 
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FIG. 2: The geodesies in equilibrium which are relevant 
to the transport of Hawking radiation from the horizon to 
the boundary for a specified \q\/ui [8, 9]. (The horizon is at 
r — nT = 1 in the coordinates of eq. (7).) Motivated by Fig. 1, 
the grey band indicates where the non-equilibrium geometry 
can not be reasonably modelled by equilibrium. 

are well known from computations of jet-energy loss in 
strongly coupled plasmas [8, 9]) propagate up from the 
horizon and skim just below the boundary for a para- 
metrically long time compared to the final temperature, 
t ~ (wcr) 1 / 4 /7rTf . (This behavior is shown in Fig. 2 for the 
equilibrium geometry.) Thus, lightlike boundary correla- 
tors at late times, t < (ujct) 1 / 4 /nTf , reflect the emission 
of non-equilibrium Hawking radiation at earlier times 
t < l/nTf. The delayed equilibration seen in Fig. 1(d) 
is consistent with a (wc) 1 / 4 expectation. 

Finally, we note two marked differences between the 
evolution of quantum fluctuations at weak and strong 
coupling. First, note that the back-reaction between the 
fluctuations and the metric was neglected in this study. 
At weak coupling this approximation is of limited use, 
since quantum fluctuations seed instabilities which grow 
exponentially, and invalidate the approximation scheme 
before the plasma thermalizes [11-15]. In the gravita- 
tional theory, however, the timescale for growth is the 
same as the equilibration timescale of the background 
metric. Thus, in the SYM plasma a nascent instability 
can grow by at most one e— folding before being damped 
by the equilibrium black hole metric. Second, note that 
the fluctuations equilibrate even in the collisionless (or 
linearized) approximation that we are using. In weakly 
coupled (but strong field) plasmas the analogous Hartree 
approximations fail satisfy the FDT at late times [15, 16]. 
These marked differences with weakly coupled plasmas 
provide a useful foil to perturbative analyses of the thcr- 
malization in heavy-ion collisions. 

In summary, we have computed the Hawking flux from 
non-equilibrium black holes. From a gravitational per- 
spective, the computation is both conceptually and tech- 
nically challenging. The physics of this novel emission 
process determines the qualitative features of the ther- 
malization of fluctuations in strongly coupled plasmas. 
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